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ABSTRACT 

This paper reports high-precision Stokes V spectra of HD 191612 acquired using the ES- 
PaDOnS spectropolarimeter at the Canada-France-Hawaii Telescope, in the context of the 
[ Magnetism in Massive stars (MiMeS) Project. Using measurements of the equivalent width 

■ of the Ha line and radial velocities of various metallic lines, we have updated both the spec- 

troscopic and orbital ephemerides of this star. We confirm the presence of a strong magnetic 
field in the photosphere of HD 191612, and detect its variability. We establish that the lon- 
gitudinal field varies in a manner consistent with the spectroscopic period of 537.6 d, in an 
approximately sinusoidal fashion. The phases of minimum and maximum longitudinal field 
are respectively coincident with the phases of maximum and minimum Ha equivalent width 
and Hp magnitude. This demonstrates a firm connection between the magnetic field and the 
k>( _ processes responsible for the line and continuum variability. Interpreting the variation of the 

• longitudinal magnetic field within the context of the dipole oblique rotator model, and adopt- 

ing an inclination i = 30° obtained assuming alignment of the orbital and rotational angular 
momenta, we obtain a best-fit surface magnetic field model with obliquity /3 — 67 ± 5° and 
polar strength = 2450 + 400 G . The inferred magnetic field strength implies an equatorial 
wind magnetic confinement parameter 77* =* 50, supporting a picture in which the Ha emission 
and photometric variability have their origin in an oblique, rigidly rotating magnetospheric 
structure resulting from a magnetically channeled wind. This interpretation is supported by 
our successful Monte Carlo radiative transfer modeling of the photometric variation, which 
assumes the enhanced plasma densities in the magnetic equatorial plane above the star im- 
plied by such a picture, according to a geometry that is consistent with that derived from 
the magnetic field. Predictions of the continuum linear polarisation resulting from Thomp- 
son scattering from the magnetospheric material indicate that the Stokes Q and U variations 
are highly sensitive to the magnetospheric geometry, and that expected amplitudes are in the 
range of current instrumentation. 
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1 INTRODUCTION 

The classification Of?p was first introduced by Walborn (1972) to 
describe spectra of early O-type stars exhibiting the presence of 
C m .44650 emission with a strength comparable to the neighbour- 
ing N in lines. Well-studied Of?p stars are now known to exhibit 
recurrent, and apparently periodic, spectral variations (in Balmer, 
He i, C m and Si in lines), narrow P Cygni or emission components 
in the Balmer lines and He i lines, and UV wind lines weaker than 
those of typical Of supergiants (see Naze et al. 2010 and references 
therein). 

The subject of this paper, HD 191612, is perhaps the best- 
studied of the 5 known Galactic Of?p stars (Walborn et al. 2010). 
Walborn (1973) was the first to note the distinctive peculiarities in 
its spectrum. Large, recurrent spectral variations were discovered 
by Walborn et al. (2003). Subsequently, Walborn et al. (2004) found 
that these variations were strictly correlated with the low amplitude 
~ 540 d variations of the Hipparcos lightcurve (Koen & Eyer 2002), 
suggesting an underlying "clock" that was proposed by Walborn et 
al. (2004) to be a binary orbit. However, the detection of a magnetic 
field in HD 191612 by Donati et al. (2006, B e = -220 ± 38 G) led 
to speculation that the observed variability was in fact rotational 
modulation, driven by the magnetic field and its interaction with 
the stellar wind. In this scenario, the magnetic, photometric and 
spectral variations may be interpretable within the context of the 
Oblique Rotator Model (ORM; Stibbs 1950). In the ORM, a large- 
scale magnetic field (typically a dipole) is "frozen" into the stellar 
plasma, and tilted relative to the stellar rotation axis. As the star 
rotates, observable quantities (e.g. the line-of-sight component of 
the magnetic field, stellar brightness, emission lines) are modulated 
according to the rotational period. 

Systematic spectral and photometric studies in the optical and 
X-ray domains have been carried out by Howarth et al. (2007) and 
Naze et al. (2007). Naze et al. (2007) analysed phase-resolved X- 
ray and optical observations of HD 191612. They find that the star 
is overluminous in X-rays relative to the canonical log L x - log L bo i 
relation (by about a factor of 5), that it appears brighter in X-rays 
when the optical emission lines are strongest, that the X-ray lines 
are quite broad and that the spectrum is dominated by a cool, ther- 
mal component. Based on those and new observations, Naze et al. 
(2010) confirmed that the X-ray emission is modulated according 
to the ~ 540 d period. These authors also performed a temporal 
variance spectrum (TVS) analysis on their optical spectra, confirm- 
ing that the most significant spectral variability was associated with 
Balmer and He i lines (Walborn et al. 2003). They noted that the 
TVS profiles of Balmer and He i lines are roughly Gaussian, of 
constant phase, and centred at the rest wavelength. The line profiles 
consist of an apparently constant absorption component with super- 
imposed variable emission extending from approximately -200 to 
200 km.s~' . In contrast to the Balmer and He lines, absorption 
lines from He n and absorption/emission lines of metals exhibit 
double-peaked TVS profiles which are more typical of binary sig- 
natures. 

From an extensive spectroscopic dataset spanning 17 y, 
Howarth et al. (2007) demonstrate unambiguously that the equiv- 
alent widths (EWs) of all variable spectral lines can be reasonably 
phased according to a unique period of 537.6 ± 0.4 d. In partic- 
ular, the striking Ha EW variation exhibits strict periodicity and 
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Table 1. Summary of stellar, wind, magnetic and magnetospheric properties 
of HD 191612. 



Spectral type 


UD1 :p - Uoip 


waiDorn ei ai. \avl\j) 


* eff (A,) 


nnn -t- i nnn 


riowann et di. ^zuu / ) 




3.5 ± 0.1 


Howarth et al. (2007) 


R* (Rg) 


14.5 


Howarth et al. (2007) 


vsin/ (kms -1 ) 


<60 


Howarth et al. (2007) 


log(L*/L ) 


5.4 


Howarth et al. (2007) 


M* (M ) 


~ 30 


Howarth et al. (2007) 


logMV7(M yr-') 


-5.8 


Howarth et al. (2007) 


Vco (fans -1 ) 


2700 


Howarth et al. (2007) 


Bi (G) 


2450 ± 400 


This paper 


BO 


67 ±5 


This paper 


t]« 


50 


This paper 


W 


2 x ur 3 


This paper 


f spin 


0.33 Myr 


This paper 



is characterised by a single, relatively sharp emission maximum 
and a broader, flatter minimum. The Ha extrema are separated in 
phase by exactly one-half cycle. When the Hipparcos photometry 
(ESA, 1997) of HD 191612 are phased with the Ha ephemeris, they 
also exhibit a coherent sinusoidal variation, with extrema located 
at the phases of the Ha extrema. Howarth et al. also established 
that HD 191612 is an eccentric double-lined spectroscopic binary 
(SB2) with an orbital period of 1542 ± 14 d and a mass ratio of 
0.483 + 0.044. Those authors examined the relationship between 
the 537.6 d and 1542 d periods, concluding that they are indepen- 
dent, and that the binary orbit has no important role in the spectral 
variability in the optical. This is also the case in the X-ray domain 
(Naze et al. 2010). In Table 1 we summarise the physical and wind 
properties of HD 191612 as reported by Howarth et al. (2007) and 
as derived in this study. 

In this paper we report systematic monitoring of the magnetic 
field and spectrum of HD 191612 within the context of the Mag- 
netism in Massive Stars (MiMeS) Project (e.g. Wade et al. 2011). 
In Sect. 3 we describe the observations obtained, and update the 
ephemerides corresponding to the 537.6 d period and the 1542 d 
period. In Sect. 4 we analyse the magnetic data using the Least 
Squares Deconvolution method, and diagnose the magnetic field. In 
Sect. 5 we examine the period content of the (variable) longitudinal 
magnetic field, and demonstrate that the data can be modeled as a 
periodic, sinusoidal signal, with a period equal to the 537.6 d pe- 
riod inferred from spectroscopy. In Sects. 6 and 7 we constrain the 
strength and geometry of the magnetic field dipole component, and 
find that the geometry constraint derived from the magnetic data is 
in excellent agreement with the geometry obtained by Howarth et 
al. (2007) from modeling of the Ha variation. In Sect. 8 we evaluate 
the capability of the light variation to constrain the stellar geome- 
try, and ultimately adopt a reference geometry in Sect. 9. Finally, in 
Sect. 10 we discuss the implications of these results for our under- 
standing of the behaviour of HD 191612, and for this general class 
of magnetic O-type stars. 



2 OBSERVATIONS 

High resolution spectropolarimetric (Stokes / and V) observations 
of HD 191612 were collected with ESPaDOnS at the Canada- 
France-Hawaii Telescope. Altogether, 24 observations were ob- 
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Table 3. Parameters of the updated Ha ephemeris obtained from combining 
the new spectroscopic data with those reported by Howarth et al. (2007). 
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tained over a time span of 1590 days. The majority of the observa- 
tions were obtained between July 2008 and October 2010; one se- 
quence was also obtained in June 2006. Typically, two observations 
were acquired per night, during 2 nights each observing semester. 
The observing cadence was derived under the tentative assumption 
that the 537.6 d spectroscopic and photometric period represented 
the rotational period of the star, and would therefore represent the 
period of variation of the magnetic field. Each spectropolarimetric 
sequence consisted of four individual subexposures taken in differ- 
ent polarimeter configurations. 

From each set of four subexposures we derived Stokes / and 
Stokes V spectra following the double-ratio procedure described by 
Donati et al. (1997), ensuring in particular that all spurious signa- 
tures are removed at first order. Null polarization spectra (labeled 
AO were calculated by combining the four subexposures in such a 
way that polarization cancels out, allowing us to verify that no spu- 
rious signals are present in the data (see Donati et al. 1997 for more 
details on the definition of AO. All frames were processed using the 
CFHT's Upena pipeline, feeding the automated reduction package 
Libre ESpRIT (Donati et al. 1997). The peak signal-to-noise ratios 
(SNRs) per 2.6 km.s~' velocity bin in the reduced spectra range 
from 400-750, with the variation due principally to weather condi- 
tions. 

The log of CFHT observations is presented in Table 2. 

3 EPHEME RIDES 

3.1 Ha 

As described in detail by Howarth et al. (2007), line profiles of HD 
191612 are variable in both EW and radial velocity (see in particu- 
lar their Figs. 1, 3 and 6). We have measured the equivalent width 
of the Ho- profile from the mean spectra acquired on each night. For 
completeness, we combined these data with the results of Howarth 
et al. (2007) to give the updated ephemeris summarised in Table 3. 
It is sensibly indistinguishable from the equivalent Howarth et al. 
result (P = 537.6 d), which we adopt throughout the remainder of 
this paper. In particular, there is no phase drift with respect to the 
previous result. 

Differences between these values and those reported in Table 1 
of Howarth et al. (2007) are negligible, and do not justify formally 
updating the ephemeris at this time. 

3.2 Orbit 

The new CFHT spectra are of very high quality, and now span a 
periastron passage in the long-period, low-amplitude spectroscopic 
orbit described by Howarth et al. (2007). Velocity measurements 
have been made from the mean CFHT spectra acquired on each 
night, and additionally from a spectrum obtained in service mode 



Table 4. Updated orbital solution. 
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on JD 2455264.76, close to the expected time of periastron passage, 
with the FIES spectrograph on the 2.56m Nordic Optical Telescope 
(388-598 nm,K = 46,000). 

The orbital solution uses the C iv doublet, /U5800, to estab- 
lish the primary's y velocity; /1/15800 plus a complex of 8 emission 
lines (N n, Si iv, and unidentified) around AA67QQ to determine all 
other orbital parameters for the primary; and a set of 6 O n lines 
(/1/14300-4700, with separations constrained to laboratory values) 
to determine K 2 (and y 2 ; see Howarth et al. 2007 for further details). 
The typical external errors on the new velocity measurements are 
~1 km.s-' for C iv and -2-3 km.s-' for AA6100 and O n. 

The revised orbit, which affords a modest improvement over 
the previous solution because of the availability of superior data at 
critical phases, is summarised in Table 4 and illustrated in Fig. 1. 

The new CFHT RV measurements are reported in Table 2. The 
FIES spectrum yields the following velocities: +17.2 km.s -1 for 
,15800, -37.9 km.s- 1 for On. 



4 MAGNETIC DIAGNOSIS 

An examination of the reduced spectra reveals weak Zeeman sig- 
natures visible in profiles of strong absorption lines during phases 
of peak Ha emission. However, the SNR of these signatures is not 
sufficient for detailed modeling. Instead, we will combine the sig- 
nal from a handful of lines in the spectrum to improve our ability 
to diagnose the magnetic field. 

Least-Squares Deconvolution (LSD, Donati et al. 1997) was 
applied to all observations. In their detection of the magnetic 
field of HD 191612, Donati et al. (2006) developed and applied a 
manually-constructed line mask containing 12 lines. In developing 
this mask, they excluded lines in emission, strong Balmer lines, and 
other features showing P Cygni profiles. To begin, we used this line 
list to extract, for all collected spectra, mean circular polarization 
(LSD Stokes V), mean polarization check (LSD AO and mean un- 
polarized (LSD Stokes /) profiles. All LSD profiles were produced 
on a spectral grid with a velocity bin of 36 km.s -1 , providing rea- 
sonable sampling of the observed mean profile and maximising the 
S/N per LSD pixel. The LSD Stokes / profile shows an extension 
to high (~ 300 km.s- 1 ) velocities in the blue wing. This asymme- 
try is similar to that observed in the LSD profiles of the Of?p star 
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Table 2. Log of ESPaDOnS observations of HD 191612. All exposures were 4800 seconds duration except those corresponding to odometer #s 851550- 
851562, which were 2240 s. Listed are the unique CFHT odometer number of the first spectrum of each 4-subexposure sequence, the heliocentric Julian 
date of the midpoint of the observation, the peak signal-to-noise ratio per 2.6 km.s~' velocity bin, the epoch and phase of the observation (according to 
the spectroscopic ephemeris of Howarth et al. 2007, Eq. 1), the evaluation of the detection level of a Stokes V Zeeman signature (DD=definite detection, 
MD=marginal detection, ND=no detection), and the derived longitudinal field and longitudinal field detection significance z from both V and N. In no case 
is any marginal or definite detection obtained in the N profiles. Also included are the Ha equivalent widths and radial velocities measured from the averaged 
spectra acquired on each observing night, and described in Sect. 3. The radial velocities and EW corresponding to odometer #s 85 1550-85 1562 were published 
previously by Howarth et al. (2007), and we include those values here for completeness. 
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HD 148937 (Wade et al. 2011) although not quite as strong. It is 
likely reflective of the inclusion of asymmetric He i lines in the 
mask. An additional dip in the blue wing is attributed to the dif- 
fuse interstellar band (DIB) located at ~ 579.8 nm, which was not 
included in the line mask but which is blended with C iv /15801, 
which is included. 

Using the^ 2 signal detection criteria described by Donati et al. 
(1997), we evaluated the significance of the signal in both Stokes 
V and in N. In no case is any signal detected in A', while signal 
in V is detected marginally (false alarm probability fap < 1(T 3 ) or 
definitely (fap < 10~ 5 ) in 15 of our 24 observations. Finally, from 
each set of LSD profiles we measured the mean longitudinal mag- 
netic field in both V and N using the first moment method (Rees 
& Semel 1979) as expressed by Eq. (1) of Wade et al. (2000), inte- 
grating from -280 to +180 km.s -1 . The longitudinal field measured 
from Stokes V is detected significantly (i.e. |z| = \B f \/cr > 3) in 11 
of our observations. In no case is the longitudinal field significantly 
detected in N. 

The results of the analysis described above are summarised in 
Table 2. 

To test the sensitivity of these results to the detailed mask 
composition, we extracted LSD profiles for two additional lines 
masks. We began using a generic line mask based on a 40000 K ex- 
tract stellar request from the VALD database, using a line depth 



threshold of 0. 1 . We sought to construct a mask that yields a mean 
LSD Stokes / line profile that exhibited as little variability as possi- 
ble, and that was as symmetric as possible, while still maximizing 
Stokes V signal. At each step of the mask development, we visu- 
alised the agreement of the LSD model (i.e. the convolutions of the 
Stokes / and V LSD profiles with the line mask) with the reduced 
spectrum, and evaluated the symmetry and variability of the LSD 
profile. We quickly realised that many lines, even if in absorption 
during some phases of the 537.6 d cycle, needed to be discarded 
from the mask due to emission at other cycle phases. After removal 
of the most significantly affected lines, we were left with a mask 
containing 26 lines, for which we adjusted the mask depths to best 
match the mean observed line depths. While none of the lines re- 
maining in the mask show obvious emission, many are asymmetric 
and vary during the course of the 537.6 d cycle. This mask is dom- 
inated by lines of neutral He, but also contains lines of ionised He, 
C, N O and Si. LSD profiles extracted using this mask - which con- 
tains more than twice the number of lines in the Donati et al. mask 
- yields LSD profiles nearly identical to those obtained from the 
Donati et al. mask. This is likely a consequence of the inclusion 
of many of the same strong He lines in both masks. Measurements 
of the longitudinal fields extracted from these LSD profiles are in 
almost perfect agreement with those obtained above. 

As a second test, we extracted LSD profiles using a more re- 
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stricted mask containing only a half-dozen metallic lines (lines of 
C, and Si). This mask yields profiles whose shapes differ sub- 
stantially from those obtained from the first two masks: they are 
narrower, and they are nearly symmetric. This is expected due to 
the exclusion of the much broader, more asymmetric He i lines. 
Notwithstanding the difference in shape, however, this mask yields 
longitudinal fields that are in good agreement (although slightly 
noisier due to the weaker and more limited sample of lines) with 
those obtained from the masks containing He lines. Similar agree- 
ment is obtained whether we employ an integration range adapted 
to the (narrower) width of the metallic lines, or if we use the (larger) 
width adapted to the masks containing He lines. 

Based on this analysis, we conclude that while the detailed 
shapes of the LSD profiles depend on the composition of the mask, 
the inferred longitudinal field of HD 191612 is relatively insensitive 
to the mask composition. We therefore based our following analysis 
on the LSD profiles extracted using the 12 line mask of Donati et 
al. (2006). 

As mentioned above, we typically obtained two observations 
of HD 191612 during each observing night (or, in one case, within 
several nights). We confirmed that no significant variability had oc- 
curred between any of these near-simultaneous observations. This 
is consistent with the known long variability period of this star. To 
improve the precision of our measurements, in our following anal- 
ysis of the magnetic variability and geometry of HD 191612 we 
have coadded and remeasured LSD profiles obtained close in time 
(i.e. on the same night or, in the case of observations 1204158 and 
1204177, within 3 nights). Ultimately, this has yielded 12 inde- 
pendent measurements of the magnetic field. These measurements, 



which will be used for the remainder of the analysis, are reported in 
Table 5. Examples of the coadded LSD profiles are shown in Fig. 
2. 

In our following analysis, we also include the observa- 
tion of Donati et al. (2006) of -220 + 38 G obtained on HJD 
2453546.55279. 



5 MAGNETIC PERIOD 

The longitudinal magnetic field measurements of HD 191612 are 
strongly inconsistent with the null field hypothesis (reduced x 2 of 
19.4). This confirms the magnetic field detection reported by Do- 
nati et al. (2006). Nor are they consistent with a nonzero but non- 
variable field (reduced^ 2 of 8.6). The longitudinal magnetic field 
of HD 191612 is therefore detected to vary on the timescales of our 
observations, with a maximum measured value of 67 ± 73 G and a 
minimum value of -514 ± 81 G. 

Including the measurement of Donati et al. (2006), the longitu- 
dinal field measurements span a total time of 1940 days. The Lomb- 
Scargle periodogram of the magnetic measurements between 1 and 
1000 days shows a single minimum with reduced^ 2 below 2.5. The 
minimum is double-peaked, with dominant peak (with reduced x~ 
of 0.5) at 539 d and a secondary peak (with reduced x 2 of 1.5) at 
623 d, both ±20 d approximately. The 539 d period is consistent 
with the well-determined spectroscopic period of 537.6 ± 0.4 d, 
while the 623 d period may be a sidelobe. Although the magnetic 
data set only spans 3.6 cycles of this length, the lack of other signifi- 
cant power in the periodogram supports a basic tenet of the oblique 
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Figure 2. LSD profiles obtained using the Donati et al. mask. Left: Stokes V, N and / LSD profiles at phase 0.979 (Br = -514 ± 81). Right: Stokes V, N and / 
LSD profiles at phase 0.601 (Be = -65 ± 63). The profiles have been scaled according to the SNR-weighted mean wavelength, lande factor and line depth of 
the mask. Note the significant variation of the strength of the Stokes / profile. 



rotator model (ORM): that rotational modulation determines the 
variability of magnetic and other observables, and that their peri- 
ods should therefore generally be the same. Indeed, when phased 
with the 537.6 day spectroscopic period (Howarth et al. 2007) the 
magnetic data display a smooth, approximately sinusoidal varia- 
tion, with a very low reduced ;f 2 of 0.5 (which may suggest that our 
error bars are slightly overestimated). We note that the data do not 
phase satisfactorily according to the SB2 orbital period of 1542 d 
(or other periods in this vicinity), yielding a best reduced^ 2 ~ 6.5). 

We therefore conclude that the variable magnetic field of HD 
191612 is best interpreted within the context of the ORM. In this 
case, the common magnetic/spectroscopic/photometric period of 
537.6 d is the rotational period of the star. The magnetic obser- 
vations therefore confirm the (until now) tentatively assumed ORM 
framework (e.g. Donati et al. 2006, Howarth et al. 2007) for inter- 
pretation of the observed variability, and the extreme slow rotation 
of HD 191612. Moreover, extending these conclusions to the cases 
of other Of?p stars in which magnetic fields have been detected 
(HD 108 [Martins et al. 2010], HD 148937 [Wade et al. 2011]) 
supports the view that all magnetic Of?p stars host strong, organ- 
ised magnetic fields, that they may well be oblique rotators, and 
that their diverse spectroscopic periods (from 7 d for HD 148937 
to perhaps 55 y for HD 108) likely represent their real rotational 
periods. 



acquisition of those observations. The phased Ha data (both new 
and old) are illustrated in the middle frame of Fig. 3. 

What is immediately apparent in Fig. 3 is that the extrema of 
the Ha equivalent width variation, as well as the Hipparcos photo- 
metric variation (lower frame), occur simultaneously with the ex- 
trema of the magnetic variation. This further strengthens the view 
implicit to the oblique rotator model that there exists a causal rela- 
tionship between the magnetic field variation, and the spectral and 
photometric variations. 

We note that Hubrig et al. (2010) report a single SOFIN obser- 
vation of the longitudinal field of HD 191612, B c = +450 ± 153 G, 
acquired on September 11, 2008, and corresponding to rotational 
phase 0.43. This strong positive longitudinal field does not agree 
very well with our observed field variation, which corresponds to 
(an essentially) consistently negative longitudinal field. 

Hubrig et al. also report the presence of very strong (up to 
15%) Stokes V signatures in the Ca n H and Na i D lines that have 
the same (positive) sign as the stellar longitudinal field they infer. 
However, these lines are interstellar features in the spectrum of HD 
191612, and consequently a linkage with the stellar field is unlikely. 
We observe no similar Stokes signatures in our ESPaDOnS spectra 
of these lines. The strong interstellar signatures might arise from 
the non-standard normalization of the Stokes V spectrum, which 
Hubrig et al. indicate to be V/I rather than the more usual V/I c ); 
see their Fig. 1. 



6 PHASE VARIATION OF THE LONGITUDINAL FIELD 

To model the magnetic field of HD 191612 within the framework 
of the oblique rotator model, we begin by phasing all of the data 
according to the Ha ephemeris of Howarth et al. (2007): 

JD =2453415.2(5) ± 537.6(4) ■ E. (1) 

The phased longitudinal field measurements are illustrated in 
the upper frame of Fig. 3. 

To verify the phasing of the magnetic data relative to the spec- 
troscopic ephemeris, we have also measured the equivalent width 
of Ha from each of the co-added spectra, using the same methods 
described by Howarth et al. (2007). The new Ha measurements 
phase perfectly with the extensive dataset published by Howarth et 
al., indicating that no significant phase drift has occurred since the 



7 STELLAR GEOMETRY IMPLIED BY THE 
MAGNETIC VARIATION 

The least-squares sinusoidal fit B [((/>) = So + 5i cos 2n((j> + tfn) to 
the phased magnetic data yields coefficients Bo = -236 ± 10 G, 
S[ = 324 ± 20 G, and 0! = 0.48 ± 0.04 corresponding to extrema of 

= +88 ± 22 G, B™ n = -560 ± 22 G. The reduced x 2 of this fit 
is 0.5, suggesting an overestimation of the longitudinal field error 
bars by a factor of ~ 1 .4, assuming a target reduced x 2 of unity. 
This may be a consequence of the relatively small number of lines 
used in the LSD line mask. 

We model the magnetic field of HD 191612 assuming an 
oblique rotating dipole, characterised by four parameters: the phase 
of closest approach of the magnetic pole (0i - 1 = -0.02, obtained 
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Figure 3. Longitudinal magnetic field, Ha EW and Hipparcos magnitude, all phased according to the 537. 6d period. The measurement of Donati et al. (2006) 
appears at phase 0.23. Solid lines are least-squares fits to the data. 



from the least-squares fit), the stellar rotation axis inclination (i) 
to the observer's line-of-sight, the magnetic axis obliquity (jS) and 
the dipole polar strength B A . The longitudinal field variation pro- 
vides constraint on two of these unknowns, typically /3 and 6 d . Fre- 
quently, the inclination i is inferred using the measured projected 
rotational velocity and period, and the inferred stellar radius. How- 
ever, in the case of HD 191612, only an upper limit on vsini is 
known (< 60 km.s~' ), and this value is clearly highly in excess of 
the true v sin i if P mt = 537.6 d. We therefore cannot infer i using 
the usual methods. 

To begin, assuming the inclination angle to be unconstrained, 
we have modeled the longitudinal field variation for i ranging 
from 20 - 80°, obtaining best fit values 13° < ji < 75° and 
2000 G< |B d | < 5000 G. For inclination angles approaching 0/90°, 
the obliquity approaches zero and the dipole strength diverges. On 
the other hand, models with i = 0° and i = 90° are not able to 
reproduce the longitudinal field variation. The general family of 
solutions is characterised by i +J3 = 95 ± 10°. 

Howarth et al. (2007) attempted to reproduce the form of the 
Ho- EW variation using two simple models: first, a centred, tilted, 
geometrically thin but optically thick Ha-emitting disc, and sec- 
ondly, a single surface spot. Both are equally capable of reproduc- 
ing the Ho- variation, and both models yield the basic geometrical 
constraint i+a ^ 100° , where a is the obliquity of the disc or the co- 



latitude of the spoQ. This constraint is consistent with that derived 
from the longitudinal field variation, substituting ft for a, i.e. that 
the a angle employed by Howarth et al. (2007) is in fact equal to 
the magnetic obliquity. In the particular case of the thin disc model, 
this implies that the disc lies in the plane of the magnetic equator, as 
would be expected in the framework of a magnetically-channeled 
wind (e.g. Babel & Montmerle 1997, ud Doula & Owocki 2002). 



8 MONTE CARLO RADIATIVE TRANSFER 
SIMULATION OF THE LIGHT VARIATION 

To leverage the diagnostic potential offered by the Hipparcos light 
curve (Fig. 3), we use a Monte-Carlo radiative transfer (RT) code 
developed by RHDT for simulating light scattering in circumstellar 
envelopes. The RT code follows much the same procedure as other 
recent codes (see, e.g., Wood & Reynolds 1999). Photon packets 
are launched from a central star and allowed to propagate through 
an arbitrary distribution of circumstellar matter (described by a 
Cartesian density grid), until they are scattered by free electrons. 
We neglect the possibility of packet absorption, since at tempera- 
tures similar to the r cff of HD 191612 the optical (Paschen contin- 
uum) bound-free cross section is negligible compared to the Thom- 



105° was stated by Howarth et al., but a revisit of the modeling by IH 
yields the slightly lower 100°. 
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Figure 4. Photometric (left) and linear polarisation (centre and right) phase variations of HD 191612 computed for various rotational axis inclinations. The 
open circles in the left-hand panel represent the Hipparcos photometry. The Hipparcos photometry of HD 191612 constrains i > 30°. 



Table 5. Averaged data used for modeling. 
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2455486.796 


3.853 


0.853 


DD 
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son cross section. Upon scattering, a ray is peeled off from the 
packet toward a virtual observer, who records the packet's Stokes 
parameters appropriately attenuated by any intervening material 
(see Yusuf-Zadeh, Morris & White 1984, for a discussion of this 
peel-off technique). A new propagation direction is then chosen 
based on the dipole phase function (Chandrasekhar 1960), and the 
packet's Stokes parameters are updated to reflect the linear polar- 
ization introduced by the scattering process. The propagation is 
then resumed until, after possible further scatterings, the packet 
eventually escapes from the system or is reabsorbed by the star. 

The initial emission direction of photon packets, relative to the 
local stellar-surface normal, is determined randomly in accordance 
with the limb-darkening law tabulated by Chandrasekhar (1960) for 
plane-parallel electron scattering atmospheres. The corresponding 
initial Stokes parameters are determined from the same tabulation. 
To improve the efficiency of the code, all packets are forced to un- 
dergo at least one scattering, using the formalism described by Witt 
(1977). 

In applying the RT code to HD 191612, we establish the cir- 
cumstellar density distribution from a 2-D (axisymmetric) MHD 
simulation of the star's magnetically channeled wind using the 
ZEUS 3-D code (see ud-Doula & Owocki 2002 and Owocki & ud- 
Doula 2004, for a general overview of our simulation approach) 



using the full energy equation as in Gagne et al. (2005), and the pa- 
rameters summarised in Table 1 . To smooth out temporal variations 
in the density, we take an average over 1,000 Ms of simulation time 
(corresponding to ~ 30 wind flow times), beginning at 300 ks after 
the initial state. Because the star rotates so slowly, the dynamical ef- 
fects of the rotation are negligible, and we can use the same density 
distribution (albeit suitably realigned to the magnetic axis) for any 
choice of the obliquity p. Electron number densities are calculated 
from the mass density under the assumption of solar composition 
and complete ionization. 

Fig. 4 shows the results from RT simulations for inclina- 
tions i = 10°, 30°, 50° and 70°, and corresponding obliquities 
P = 85°, 65°, 45° and 25° (these are the values mandated by the 
i + P = 95° relation; see Sec. 7). The left-hand panel shows the 
light curves over one rotation cycle, for the four choices of incli- 
nation. All models correctly reproduce the distinctive shape of the 
Hipparcos light curve. However, only the ( > 30° models are able 
also to match the observed ~ 4 mmag amplitude, allowing us to rule 
out the ( = 10° case. (We stress that the amplitudes of our model 
light curves are not free parameters, but rather predictions of the 
MHD and RT simulations for the adopted stellar properties). 

From the light curves alone, it is not possible to distinguish 
between the i = 30°, 50° and 70° models. However, the centre and 
right-hand panels of Fig. 4, which plot the models' Stokes Q and 
U parameters over one rotation cycle, clearly reveal that the linear 
polarization variations of the star are sensitive to inclination. More- 
over, the presence and variation of the linear polarisation would 
implicitly confirm the presence of a flattened, disc-like structure. 
This finding establishes a strong case for obtaining fresh phase- 
resolved linear continuum polarization observations of HD 191612. 
Such observations - which require a precision of order 0.01% - are 
within range of current instrumentation (e.g. Carciofi et al. 2007). 



9 TENTATIVE REFERENCE GEOMETRY 

Although the inclination i is not strongly constrained by the mea- 
surements in hand, the data provide a relatively precise diagnosis 
of the magnetic parameters p and for a given value of i. Given 
our desire to present a reference model for the magnetic field of 
HD 191612, we proceed as follows. Hypothesising that the compo- 
nent masses in the binary system are more or less normal for their 
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Figure 5. Reduced x 2 landscape of /3 vs Bj for fixed inclination i = 30°. 
Contours correspond to l-5<x limits obtained from fitting the (two free pa- 
rameter) oblique magnetic dipole model to the phased, averaged longitudi- 
nal field variation. We obtained B d = 2450 ± 400 G and = 67 ± 5°, where 
the error bars are the formal (lcr) uncertainties. 

spectral types (08fp for the primary, B1V for the secondary), this 
suggests that sin(7 OTb ) is very likely within the range 0.4-0.6 (Fig. 9 
of Howarth et al. 2007). Speculating that the orbital and rotational 
angular momenta are aligned, we tentatively adopt i = i orb ^ 30° 
as a first guess at the rotation axis inclination. We note that this 
inclination is consistent with the observed Hipparcos photometric 
variation according to the modeling described in the previous sec- 
tion. 

Modeling the longitudinal field variation using a numerical 
implementation of the dipole oblique rotator model, for i = 30°, 
we obtain B A = 2450 ± 400 G and p = 67 + 5°. For illustration, 
the reduced x 2 landscape for p vs. B A corresponding to i = 30° is 
presented in Fig. 5. 

The results of our magnetic field modeling are in remarkably 
good agreement with the strength and geometry proposed by Do- 
nati et al. (2006) based on their single measurement of the lon- 
gitudinal field, and guidance provided by the Ha equivalent width 
variation. In particular, they proposed that the geometry of the mag- 
netic field is constrained according to i + p = 90° (essentially the 
constraint we derive here from a complete modeling of the phase 
variation of the longitudinal field). Their inferred dipole strength 
of 1.5 kG derived for i = p = 45° is very similar to that which we 
derive for this same geometry, although somewhat weaker than that 
derived for our i = 30° reference geometry. 



10 DISCUSSION 

Our new magnetic observations of HD 191612 confirm the exis- 
tence of an organised magnetic field in the photosphere of this 
Of?p star. We have demonstrated that the longitudinal magnetic 
field varies significantly, with a period that is in agreement with 
the spectroscopic period of 537.6 d. The observed phenomenology 
is consistent with an oblique magnetic rotator, for which the spec- 
troscopic period is the stellar rotation period. 

The Ho- equivalent width and Hipparcos photometric measure- 
ments vary in phase with the longitudinal field. This strongly sug- 
gests a causal connection between their variability and the longitu- 
dinal field. In particular, Ha emission from, and electron scatter- 
ing in, a disc of plasma located in the magnetic equator are able to 



qualitatively and, to some extent quantitatively, reproduce all ob- 
servables. 

Modeling the longitudinal field variation and adopting a rota- 
tion axis inclination to the line of sight i = 30° (which assumes 
alignment of the spin and orbital angular momenta of HD 191612, 
and which is consistent with the lower limit on i implied by our 
modeling of the light variations), we derive the dipole strength and 
obliquity B A = 2450 ± 400 G and p = 67 + 5° . Using the stellar and 
wind parameters reported by Howarth et al. (2007), we compute 
the wind magnetic confinement parameter rj, = B^R 2 / Mv^ = 50 
and rotation parameter W = v eq /v cr it = 2 x 10~ 3 (e.g. ud Doula et al. 
2008). This places the Alfven radius at about Rah = r] 1 , 14 = 2.2 R,. 
The Kepler (or corotation) radius is located much further away, 
at about R Ksp = W~ 2 ' 3 = 60 R,. As pointed out by ud Doula et 
al. (2008), for any material trapped on magnetic loops inside the 
Kepler radius, the outward centrifugal support is less than the in- 
ward pull of gravity; since much of this material is compressed 
into clumps that are too dense to be significantly line-driven, it 
eventually falls back to the star following complex patterns along 
the closed field loops. Hence, for HD 191612, all magnetically- 
confined plasma (i.e. all wind plasma located inside the Alfven 
radius) is unstable to this phenomenon. Therefore the modulation 
of the Ha emission (and presumably much of the other line vari- 
ability, and the photometry as well) is a consequence of a statis- 
tical overdensity of plasma near the magnetic equatorial plane in- 
side Rah. However, this dense plasma is not supported either by 
the radiative or the centrifugal forces, and thus exhibits a rela- 
tively short residence time. The "disc" of HD 191612 is therefore 
physically distinct from the rotationally-supported, rigidly-rotating 
magnetospheres of rapidly-rotating magnetic stars such as <r Ori 
E (Townsend, Owocki & Groote 2005) or HR 7355 (Oksala et al. 
2010, Rivinius et al. 2010). 

Naze et al. (2007) highlight basic differences between the X- 
ray spectrum of HD 191612 and those of the known magnetic hot 
stars t Sco and 6 l Ori C, both of which (like HD 191612) are char- 
acterised by Rah « Rkc P - In particular, the X-ray spectrum of HD 
191612 displays a relatively soft spectrum with few high-ionisation 
lines, broad (~ 2000 km.s~' ) lines, and lack of bremsstrahlung con- 
tinuum. In contrast, r Sco and Q l Ori C show narrow X-ray lines 
(always below 1000 km.s~' , and typically a few hundred km.s~' ), 
hard spectra with a dominant hot component and high ionisation 
lines. 8 l Ori C furthermore exhibits bremsstrahlung continuum, 
clearly indicating that the emission measure distribution is domi- 
nated by high temperature plasma. They proposed that HD 191612 
is an intermediate case between the magnetic OB stars and the 
"typical" O stars. The results presented here demonstrate that the 
magnetic and magnetospheric characterstics of HD 191612 (dipo- 
lar field topology, field strength of order kG, intermediate obliquity, 
j], and W parameters) are rather similar to those of 9 [ Ori C. Hence 
the differences in X-ray characteristics of these stars must now be 
understood in light of our confirmation of the magnetic oblique ro- 
tator model for HD 191612. Given the relatively high quality of the 
current X-ray and other observations of Of?p stars, this will likely 
require guidance to theory provided by a pan-spectral approach to 
a large and diverse sample of magnetic OB stars. 

The inferred rotation period of HD 191612 - approximately 



2 If we use the theoretical mass-loss recipe of Vink et al (2000), we obtain 
an unclumped mass-loss rate of log(M) = -6. 1 (twice smaller than that 
found by Howarth et al. (2007), resulting in a confinement parameter of 
f/, ~ 100 and an Alfven radius of i?Aif = 3.2/?,. 



10 



G.A. Wade et al. 



1.5 years - is extremely slow. Computing the magnetic braking 
spindown time according to Eq. (25) of ud Doula et al. (2009), 
using the parameters in Table 1, we obtain 0.33 Myr. If the age of 
HD 191612 is 3-4 Myr as suggested by the analysis of Donati et 
al. (2006), we conclude (as did Donati et al. 2006 and ud Doula et 
al. 2009) that magnetic braking alone could be responsible for in- 
creasing the rotational period of HD 191612 from that of a typical 
ZAMS early O-type star to its current value. We note however that 
this model of rotational spindown does not take into account the 
(strongly) changing stellar moment of inertia, mass loss character- 
istics and magnetic characteristics during the evolution of the star 
from the ZAMS. In particular, if the surface dipole magnetic field 
of HD 191612 is currently of order 2 kG, under the assumption of 
magnetic flux conservation it may have been larger than 30 kG at 
the ZAMS (or possibly even larger than that, given the results of 
Landstreet et al. 2008). We conclude that more sophisticated mod- 
els taking into account the time-dependence of the spindown must 
be developed to accurately compare observed magnetic stars with 
spindown predictions. 

The identification of this period with the period of stellar ro- 
tation assumes the validity of the ORM. Such a long rotational pe- 
riod implies an equatorial rotational velocity of order 1 km/s, which 
would not be detectable in the photospheric lines of this star. Sim- 
ilarly, the ~ 55 y period inferred for HD 108 (Martins et al. 2010), 
if it is in fact the rotational period, also implies a negligible v sin ('. 
The remarkable similarity in shape of some of the least variable, 
presumably most "photospheric" lines (those least modified by the 
wind or magnetosphere, e.g. O in 5590) of these stars is consistent 
with this idea. On the other hand, those same lines in the spectrum 
of HD 148937, which according to the ORM should be rotating 
much faster (with a period of 7 d, with v eq ^ 108 km.s~' ; Wade et 
al. 201 1), also show profiles that differ negligibly from HD 191612. 
This suggests that the profiles of those lines contain no significant 
information about the stellar projected rotational velocity. 

Our confirmation of the oblique rotator model for HD 191612 
provides a reasonable basis for proposing that other Of?p stars in 
which magnetic field is detected (HD 108 - Martins et al. 2010, HD 
148937 - Wade et al. 201 1) are also oblique rotators, with their vari- 
ation periods being their rotational periods. We have proposed that 
the general class of Of?p stars may be the first category of O-type 
magnetic oblique rotators. It should be a priority to investigate the 
remaining Galactic Of?p stars (NGC 1624-2 and CPD-28°2561; 
Walborn et al. 2010) to seek to detect magnetic fields, and investi- 
gate their rotation and wind confinement. 
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